FR-DE Adequacy¶

In this example, we will test hadar on a realistic (yet simplify) use case. We will perform adequacy between France and Germainy during one day.

import pandas as pd
import numpy as np
import hadar as hd

Import simplify dataset¶

fr = pd.read_csv('fr.csv')
de = pd.read_csv('de.csv')

Build study¶

study = hd.Study(['fr', 'de'], horizon=48)

France loves nuclear, so in this example most of production are nuclear. France has also a bit of solar and when needed country can turn on/off coal generator. We want to optimize adequacy by reduce CO2 production. Therefore:

solar is the cheaper at 10
then we use nuclear at 30
and coal at 100

study \
    .add_on_node('fr', data=hd.Consumption(name='load', cost=10**6, quantity=fr['cons'])) \
    .add_on_node('fr', data=hd.Production(name='solar', cost=10, quantity=fr['solar'])) \
    .add_on_node('fr', data=hd.Production(name='nuclear', cost=30, quantity=fr['nuclear'])) \
    .add_on_node('fr', data=hd.Production(name='coal', cost=100, quantity=fr['coal']));

Germainy has stopped nuclear to switch from renewable energy. So we increase solar and eolien production. When renewable energy are off, Germainy need to start coal generation to match its consumption. Like for France, we want to minimize CO2 production:

solar at 10
eolien at 15
coal at 100

study \
    .add_on_node('de', data=hd.Consumption(name='load', cost=10**6, quantity=de['cons'])) \
    .add_on_node('de', data=hd.Production(name='solar', cost=10, quantity=de['solar'])) \
    .add_on_node('de', data=hd.Production(name='eolien', cost=15, quantity=de['eolien'])) \
    .add_on_node('de', data=hd.Production(name='coal', cost=100, quantity=de['coal']));

Then both side links are set with same cost at 5. In this network, Germany will be import from nuclear french before to start coal. And France will use germain coal to avoid any loss of load.

study \
    .add_link(src='fr', dest='de', cost=5, quantity=4000) \
    .add_link(src='de', dest='fr', cost=5, quantity=4000);

optimizer = hd.LPOptimizer()
res = optimizer.solve(study)

agg = hd.ResultAnalyzer(study, res)
plot = hd.HTMLPlotting(agg=agg,
                       unit_symbol='MW', # Set unit quantity
                       time_start='2020-02-01', # Set time interval
                       time_end='2020-02-02')

plot.stack(node='fr')

plot.stack(node='de')

Hadar found a loss of load near 6h in Germainy and import from France. Then France had a loss of load, and Hadar exports to France.